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First results of an ongoing high statistics study of the colour flux distribution around static quark sources in 
SU{2) gauge theory are presented. The flux tube profiles and widths have been investigated for several quark 
separations at /J = 2.5 and (3 = 2.74. The results are tested against Michael's sum rules. 
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1. INTRODUCTION 

Recently, high statistics studies of the static qq 
potential in SU (2) and SU (3) gauge theories have 
been performed . A linear long range potential 
with the universal subleading —tt/{12R) correc- 
tion (for quark separations Ra > 0.5 fm), as pre- 
dicted by the string picture is observed. More- 
over, the gap between the ground state potential 
and the first excited state potential for large R is 
consistent with the string value of vr/i? The 
hope for a deeper understanding of the underlying 
dynamics of the confinement mechanism is one 
of the motivations for the present investigation. 
Effective models like a bosonic string model 
based on a Nielsen-Olesen szenario of magnetic 
confinement Q can be tested. 

Measurements of colour flux distributions are 
extremely difficult because the relevant correla- 
tors of pure glue quantities exhibit large statisti- 
cal fluctuations. For this reason, previous stud- 
ies (e.g. Refs. d,^) have only been qualitative. 
In view of the great progress achieved recently 
in hardware and lattice techniques (improved op- 
erators, noise reduction, updating algorithm), it 
appears to be of considerable interest to perform 
really reliable quantitative studies of field distri- 
butions. I will report on first results of an ongoing 
project with this aim. 

2. SIMULATION 

We study lattice volumes of 16'* at = 2.5 (a « 
0.085 fm) and 32^ at /3 = 2.74 (a w 0.041 fm). 



Thus, the physical volumes are (1.3 fm)'* in ei- 
ther case. The scale has been extracted from the 
string tension value ^/a = 440 MeV. For the up- 
dating a hybrid of heatbath and overrelaxation 
has been used (details in Ref. @]). Measurements 
of the relevant operators have been taken every 
100 sweeps. At the two (3 values, 8640 and 670 
such statistically independent measurements have 
been taken, respectively. For safety, the data was 
binned into blocks of five prior to the analysis. 

The chromoelectric (and -magnetic) field distri- 
butions in presence of two static quarks have been 
measured for various quark separations r = Ra. 
The quarks are generated at time t = 0, and anni- 
hilated at time t = Ta. This is realized by use of 
Wilson loops W{R,T). A variant of APE smear- 
ing ^ has been applied to the spatial parts of 
the Wilson loops in order to increase the ground 
state overlaps. All overlaps have been found to 
be larger than 95% in the potential analysis. 
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Figure 1. Comparison between smeared and un- 
smeared operators for e(0) at /? = 2.5, i? = 4. 
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The Maxwell field strength tensor is given by 
F^^^ = 2TyTI, « 2f3/a^S^, with S^, = I - 
-^ReTr {U^^) being the plaquette action. Due to 
its locality, 5^;^ undergoes a multiplicative renor- 
malization which, in the mean field approxima- 
tion, is cancelled out in the combination PS^^. 

We measure the difference between this oper- 
ator in presence of the qq pair, separated by a 
distance r, and its vacuum expectation for vari- 
ous times T: 



p(R.T)f._ {WiR,T)S,An,T/2)) 
W- {W{R,T)) 

The (squared) field components are given by 



--^<i^JkP]k (n) 



(1) 

(2) 
(3) 



for T large. From these fields, the energy, £(x), 
and action densities, cr(x), are calculated: 



e(r)(x) 

CT(^)(x) 
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- E^,)(x)±B^,)(x: 



(4) 



Since < in the chosen Minkowski notation, 
the energy densities are very small and, thus, ex- 
tremely difficult to measure. In the actual simu- 
lation, various electric/magnetic plaquettes have 
been averaged, in order to obtain an operator that 
is symmetric around the lattice site (n, r/2). 

Early investigations have failed to identify the 
asymptotic plateau (Eqs. (2,3)). This situation 
is greatly improved by smearing techniques, as 
can be seen from Fig. where a comparison 
of e(4a)(0) (at (3 = 2.5) from smeared and un- 
smeared Wilson loops is made. Also, the statis- 
tical errors are substantially reduced (note the 
logarithmic scale). The coordinate system is cho- 
sen such that the quark source resides at the 
spatial coordinate (0, 0,i?/2), and the antiquark 
at (0,0, — i?/2). The fields have been measured 
at all points n with nj, being varied along the 
whole lattice axis, while the transverse distance 
71 = {ni,n2) is taken along the directions (1,0), 
and (1, 1) up to |n_L| = 6. 



/li/f 



A j i 



Figure 2. The action density cr(n) dX (3 — 2.5, 
r = 8a « 0.7 fm. 

3. RESULTS 

In Fig. ^, the measured action density (at /3 = 
2.5, r = 8a « 0.7 fm) is displayed. The resolu- 
tion of our lattice allows us to see smooth physical 
structures without any recourse to fancy interpo- 
lations. Note that the mesh is not equidistant in 
the perpendicular direction because the off-axis 
separations n±^ cx (1,1) are included. 

It is instructive to investigate the field com- 
ponents parallel (i?j^ — El) and perpendicular 

{E\ ^ El + E^) to the flux tube. Around the 
middle of the confining string (n = 0) < ~ 
E^ ^ is observed. Nice scaling behaviour of 
the colourfields is found between the two /?- values. 
Within statistical accuracy, continuum rotational 
invariance is restored. 

The energy and action distributions in the 
central perpendicular plane [n^ — 0) have 
been investigated. Gaussian fits e{xj_,0) = 

exp{—x\/b'^) and dipole fits e{x±,Q) = 
c^l{x\ -1-6^)'^ have been performed. In addi- 
tion, we have numerically integrated the densities 
within a circle with the cut-off radius x-casa^ 



(5) 



a^max fias been varied in order to identify a 
plateau. The width is defined by p^(r) = 
Ie{x\, r)//e(l, r). For the above fit functions one 
obtains = b^. For r < 0.4 fm the dipole 
fits yield better values. For larger qq sepa- 
rations both turn out to be acceptable but 
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the Gaussian results are closer to the numerically 
integrated values. The same holds true for p„. 
The electric field 'EP'{x^_, 0) is expected to have a 
Gaussian shape within the confining string |^ . 



energy ^ 
action r 
R^2/4 



Figure 3. Widths of the energy and action flux 
tubes p^y^ (numerically integrated) as a function 
of in lattice units. 

p1 and are plotted against E? in Fig. |[ 
Both, the action, and the energy density widths 
increase with R. In the Coulomb-dominated re- 
gion of the potential (i? < 0.4 fm) a dipole 
behaviour p|,2(^) = i?^/4 is expected. This 
straight line is plotted together with the data. 
For large separations we expect a logarithmic de- 
pendence from the string picture The actual 
lattice volumes are still too small to investigate 
this effect. 




Figure 4. Comparison between the (Gaussian 
and numerically) integrated perpendicular energy 
density and the interquark force (Eq. 6). 



The integrated energy density should equal the 
potential between the quarks up to an additive 
self energy 1^ (and the multiplicative field renor- 
malization which is approximately one for the lat- 
tice field operators chosen). Under the assump- 
tion that on increasing the qq separation, the 
additional binding energy is completely localized 
near the center of the two the sources (which is 
expected as soon as an effectively one dimensional 
string has been formed), a differential form of this 
sum rule can be derived: 

F{R) = -a^Jd^x^e^na){x±,0) , (6) 

where F{R) = —dV{R)/dR is the interquark 
force (in lattice units). This relation is tested 
in Fig. ^. Agreement is found even for i? = 2, a 
distance, where the Coulomb contribution to the 
potential clearly dominates. This is reflected in 
an extreme stability of the shapes of the fleld dis- 
tributions near the sources against increasing the 
quark separation. 

Till now, we have only realized physical quark 
separations smaller than 0.7 fm, where a one di- 
mensional string is not yet really formed. Thus, 
small differences between the fields at n3 = and 
ns — 1 are still visible. Work on a lattice of linear 
extent L « 2.7 fm is in progress 0. 
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